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The vector analyzing power A„ is calculated for elastic electron scattering from a variety of spin 
zero nuclei at energies from 14 MeV to 3 GeV. Time reversal symmetry insures that An vanish in first 
Born approximation. Therefore A„ depends on Coulomb distortions and can be large for scattering 
from heavy nuclei. The vector analyzing power is a potential source of systematic error for parity 
violation experiments. We find that A„ — —0.361 ppm for the kinematics of the Parity Radius 
Experiment (PREX) involving 850 MeV electrons scattering at six degrees from '^"^Pb. This is 
comparable to the parity violating asymmetry. However for HAPPEX He involving 3 GeV electrons 
scattering on *He we find that An is very small. 

PACS numbers; 



I. INTRODUCTION 

Recent measurements of the vector-analyzing power 
An in polarized electron-proton scattering probe nucleon 
structure and test theory Q, |^ . The analyzing power 
is a parity even, time reversal odd correlation between 
the electron spin and the target and electron momenta. 
The time reversal symmetry of Quantum Electrodynam- 
ics insures that A„ vanishes in first Born approximation. 
Therefore, a nonzero measurement of An is a direct obser- 
vation of two (or more) photon exchange contributions. 
These two-photon contributions are very interesting. For 
example, the discrepancy in the ratio of the electric to 
magnetic form factors of the proton G^/G^, obtained 
via polarization transfer and Rosenbluth separatioUjinay 
be due to the neglect of two-photon contributions 0, Q . 

The analyzing power for electron-proton scattering is 
of order the fine structure constant a. In contrast, elec- 
tron scattering from a heavy nucleus has large Coulomb 
distortion effects and these give analyzing powers of order 
Za where Z is the nuclear charge. These large analyzing 
powers are used to measure the polarization of low en- 
ergy electrons |6| . At higher energies they can be a source 
of systematic errors for parity violation experiments. If 
a primarily longitudinally polarized electron beam has a 
small transverse polarization, this may couple with the 
analyzing power and generate a false asymmetry. 

The Parity Radius Experiment (PREX)0, aims 
to measure the neutron radius of ^"^Pb accurately and 
model independently via parity violating electron scatter- 
ing. However, the analyzing power for scattering from Pb 
is relatively large because of the large charge. This could 
produce a potentially serious systematic error. There- 
fore, PREX may measure the analyzing power to control 
the error. In this paper we calculate the expected A„ for 
PREX kinematics. In addition the HAPPEX He exper- 
iment is measuring parity violating electron scattering 
from "^He at 3 GeV in order to determine strange quark 
contributions to the electric form factor of the nucleon 



0. This experiment will also measure An from ^He. 

There have not been many calculations of An at higher 
energies. In this paper we accurately calculate analyzing 
powers for electron scattering from a variety of spin zero 
nuclei over a large range of laboratory energies from 14 
MeV to 3 GeV. We solve the Dirac equation to calculate 
the analyzing power to all orders in Za. This turns out to 
be a substantial numerical challenge as discussed in Sec- 
tional In between photon exchanges, the nucleus can be 
in a variety of excited states. However, only the ground 
state generates a coherent contribution proportional to 
Za. The other excited states may only contribute at or- 
der a. Therefore the major approximation of this work 
is to assume that the nucleus remains in its ground state 
at all times. This allows the scattering to be described 
with a simple Dirac equation. 



II. FORMALISM 

In this section we describe the numerical procedure 
used in the code RUNT for calculating An. First we 
discuss the treatment of recoil. Next we discuss the nu- 
merical techniques necessary to calculate the scattering 
amplitude with enough accuracy to extract An. Since An 
is often small this is a very demanding numerical calcu- 
lation. 

We wish to describe electron scattering from a spin 
zero nucleus. We include recoil or center of mass cor- 
rections using the procedure of ref. lO]. This approach 
ends up with a Dirac equation in the relative co-ordinate, 
and an effective kernel which is a recoil factor (similar 
to a reduced mass in the non-relativistic picture) times 
the potential. Corrections can be made to this picture, 
from higher order effects, but they are not included here. 
With this approach, and our approximation that the nu- 
cleus remains in its ground state, the scattering can be 
described by the Dirac equation including the Coulomb 
potential Vc{r) from the ground state charge distribution 
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of the target. 

The time-independent Dirac equation is, 

[-la . V + /3me + ^ ^* K(r) - E,]^{r) = 0, (1) 

where Et (i?e) is the total energy (kinetic + rest mass) 
of the target (electron) in the centre of momentum frame 
and TOe is the electron mass. 

Since for large values of r, Vcif) goes like 1/r our phase 
shifts come from matching our radially integrated wave- 
functions to regular and irregular Dirac-Coulomb func- 
tions in the standard way. 

To calculate the scattering observables from the phase 
shifts there are two scattering amplitudes, 

oo 

/(0)=^A,PKcos0), (2) 

i=0 



electrons that are initially polarized transverse to the 
scattering plane. 

The problem of calculating the analyzing power for 
electron scattering by numerically solving the Dirac equa- 
tion is not an easy one. The analyzing power is very small 
for light nuclei (of the order 10^^^), so calculating it in- 
volves calculating /(0) and g(ff) to a very high precision. 
For heavier nuclei the cross section drops rapidly, so cal- 
culating an analyzing power from amplitudes that have 
already dropped several orders of magnitude is almost as 
difficult. 

In this paper we have overcome the problem using sev- 
eral techniques which are modifications to the standard 
approach. We present these techniques to help anyone 
trying to reproduce our results. 

A. Radial Integration Modifications 



.g(0)-^i?,P,i(cos0), 



where, 



and 



Ai = (? + l)?7++/?7; , 



(3) 

(4) 
(5) 

(6) 



and k is the electron's momentum in the centre of mo- 
mentum frame. 

The phase shift S{l,j) contains a part that would be 
obtained from scattering from a pure Coulomb poten- 
tial(calculated analytically), plus a piece that is obtained 
numerically from matching to Dirac Coulomb wavefunc- 
tions outside the charge distribution (typically from 10 
to 18 fm from the centre of the nucleus), and has the 
property that it goes rapidly to zero with increasing I, 
whereas the pure Coulomb phase shift goes to zero so 
slowly that the series in the expressions for f{6) and g{9) 
diverge and converge painfully slowly respectively. The 
standard approach to summing this series is well known, 
being first presented in ref fl li| . 

From the amplitudes we can calculate the cross section 
and analyzing power. 



da/dn^\fi0)\^ + \gi0)\' 



and 



An — 



2Im[/*(g)g(g)] 
urn' + 19(0)1' 



(7) 



(8) 



The analyzing power gives the fractional difference in 
cross section for scattering to the left versus right for 



The code RUNT is used for most of the calculations 
T^. The results are checked with an independent code 
ELASTIC M.,JM- RUNT uses the Numerov algorithm. 
At each radial step of length h this introduces an error 
that goes like h^, and after integrating over a certain dis- 
tance, the accumulated error goes like h'^. In practice, 
using a small enough step size to get the required preci- 
sion involved very large arrays, and we also found that 
round-off (or truncation) error of the double precision 
numbers prevented us achieving the precision we needed. 

To use the Numerov algorithm we convert the coupled 
first order Dirac equations into a single second order dif- 
ferential equation. 



where 



f{x) = -k' + U{x) + 



1(1 + 1) 



(9) 



(10) 



and the Schrodinger-equivalent potential U(x) is, 
U^(x) 



U(x) = —( 



r Ee + irie - Uc(x) 



3, U^(x) ,2 



+ 2E,Uc(x) - Uc(x)' 



U':(x) 



4:^Ee+me-Uc(x)' 2(Ee + me-Uc(x))' 



(11) 



Here Uc(x) = EtVc(x)/{Et + E^), U',{x) = dUc(x)/dx, 
and K is the usual Dirac quantum number. 

The Numerov method is a three point formula: 

2 + 10T(x) , , 

w(x + h) = — —\^w(x)-w(x-h). (12) 

1 - T{x) 

where T(x) = h'f (x)/l2 and. 



w(x) = (1 - T{x))y(x). 



(13) 
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It is self-starting for I ^ 1 and requires no power series for 
the first step. However for p- waves / = 1 it does not start 
by itself. We tried several ways of starting the integration 
but they all ended up resulting in p-wave phase shifts 
less accurate than the other phase shifts. Empirically wc 
found the best way, which was fortunately good enough, 
to be as follows. 

We make use of the fact that near the origin, any phys- 
ical potential will be an even function of r and therefore 
can be approximated as constant. When U{x) is constant 
the solution is, 

y{x)^Ajiik'x) (14) 

Where 

k'"^ ^k^ -Uih). (15) 

Thus for the first step we take 

y{2h)^y{h)n{2k'h)/h{k'h) (16) 

where ji{x) is a regular spherical Bessel function. Thus 
the first step is done this way, and then wc use use Eq. 
1121 to integrate out. This helped, but still did not result 
in a method that was accurate enough, however. 

The next step was to recognize that the kinetic en- 
ergy of the electrons is much higher than the potential 
they feel. Thus we are calculating barely perturbed plane 
waves. Under such circumstances we used two of the 
techniques discussed in ref . 

The first technique is to modify Eq. El so that it be- 
comes exact in the limit of constant f (x) . This results in 
the formula, 

w{x + h) = 2 cosh(\/l2r)w(a;) - w{x - h) (17) 

This formula still has an error that goes like ft,^, but it 
is smaller than in Eq. and this allowed us to gain 
at least two more significant figures, but this still wasn't 
accurate enough. 

The second technique is to look at where Eqs. IT^ 
and / or ITTl come from, and see that the errors involved in 
using these formulae to integrate over a fixed range (typi- 
cally ^ 18 fm in our case), go like cih'^+C2h^+c^h^ + .... 
Under such circumstances Richardsonian extrapolation 
works very well. We found in practice for this problem 
that two iterations of the Richardsonian method were op- 
timal (two iterations gained us the accuracy we needed, 
and more than this did not help as we were then at the 
round-off error limit). 

With these techniques, we found that our phase shifts 
could be calculated to the desired accuracy by match- 
ing to Dirac Coulomb wave functions calculated by the 
method of Barnet and Steed 0| . 

B. Phase-Shift Summation 

The standard approach in electron scattering, due to 
Yennie, Ravenhall and Wilson jl^ is to relate the (diver- 
gent) sum over Legendre Polynomials in equations |21 and 



|31for the scattering amplitudes /(0), g{9) with a sum over 
Legendre polynomials for the functions (1 — cos(0))"/(0) 
and (1 — cos{9))^ g{9) . These new sums converge more 
quickly, but to extract the original amplitudes we need to 
divide by (1 — cos(0))" which blows up near — 0. Thus 
we have to do our sums very very accurately at forward 
angles, n is some small integer, usually taken as 2,3,4 or 
5. This method works very well away from the forward 
angles, but is problematic there. 

Since we want to calculate accurately near forward an- 
gles, we use a different mathematical device, inspired by 
the physics of screening. In theory electron scattering 
from a point Coulomb potential gives a divergent partial 
wave series, but in practice the screening of the electron 
clouds will (for a neutral atom target), cause the partial 
wave series to converge. 

We crudely model this by multiplying the partial wave 
amplitudes Ai and Bi in equations 01 [S] by a convergence 
factor exp(— (//a)^) where a is taken to be large enough 
that our answer does not change significantly if a is in- 
creased further. 

In practice we summed these term up to about I = 6a 
(by which time the convergence factor has become very 
small) , and used a value of a somewhere between 500 
to 500,000 (depending on electron energy and target 
charge). We do caution that this method is NOT as 
numerically accurate at large angles as the method of 
Yennie, Ravenhall and Wilson. The two methods com- 
plement each other, and there is a large overlap at middle 
angles where they give the same result. 



C. Tests of Calculation 

We list some of the tests of our results. First, our re- 
sults were obtained using two independent codes RUNT 
and ELASTIC. Second we reproduce the analyzing pow- 
ers and cross sections for scattering with 14 MeV elec- 
trons presented in ref. 17]. Finally, for light (enough) 
electron mass mg, An is proportional to m,,. Should 
the electon mass get lOx heavier, and the centre-of-mass 
wave-vector remain the same, then the analyzing power 
should get lOX bigger. We explored this scaling by mul- 
tiplying the electron rest-mass by 10, 100 and 1000 and 
re-adjusted the incident electron kinetic energy so as to 
obtain the same centre of mass wave vector as was ob- 
tained by having the correct electron mass. We tested 
this for Lead at 200 MeV and for He at 3GeV (ignor- 
ing recoil effects) and found that provided the incident 
kinetic energy was not adjusted by more than 20% the 
scaling of the analyzing power worked to better then 1 %. 
This test is especially re-assuring since numerically cal- 
culating analyzing powers using lOOx heavier electrons is 
much easier. All calculations presented here do use the 
physical electron mass, but others trying to reproduce 
the numbers may find it easier to use a heavier electron 
mass. 
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III. RESULTS 



In this section we present results for A„ for scattering 
from a variety of spin zero nuclei at energies from 200 
MeV to 3 GeV. The experiment SAMPLE has measured 
An for electron-proton scattering at 200 MeV and scat- 
tering angles from 130 to 170 degrees|l|- It is interesting 
to compare this e-p result to scattering from ^"*Pb. In 
Fig. 1 we show our e— ^°*Pb cross section. We use the 
experimental charge density of ref. |3| ■ Figure 2 shows 
our prediction for An- For example, at a scattering angle 
of 146 degrees we find An — —890 ppm. This is very 
large compared to the SAMPLE result of -15.4 ± 5.4 
ppm for scattering from the proton. Our result for Pb 
is much larger because of the large charge of Pb which 
greatly increases Coulomb distortions. Also the proton 
An is somewhat reduced by a magnetic moment contri- 
bution that is absent for scattering from a spin zero 
nucleus. Figure 3 shows An at angles near 180 degrees 
where An is very large « 0.2. However, the elastic cross 
section is very small. Note that at lower energies, Mott 
polarimeters use this back angle peak in An- 
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FIG. 1: (Color online) Differential cross section for elastic 
electron scattering versus scattering angle. The solid line is 
for ^"^Pb at 200 MeV while the dashed line is for "^He at 3 
GeV. 



The Parity Radius Experiment (PREX) should mea- 
sure the parity violating asymmetry A for elastic electron 
scattering from ^°*Pb at 850 MeV and scattering angles 
near six degrees 0- The goal of PREX is to determine 
the neutron radius to 1 % The weak charge of a 
neutron is much larger than that of a proton. There- 
fore A is very sensitive to the neutron density. However, 
Coulomb distortions reduce A by approximately 30 % 
and are crucial for the interpretation of the experiment 
[T3|. A measurement of the vector analyzing power An 
directly probes these distortions since An is driven by 
Coulomb distortions and it vanishes in first Born approx- 
imation. Figure 4 shows our predictions for An at 850 
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FIG. 2; The vector analyzing power A„ for elastic electron 
scattering from '^"^Pb versus scattering angle at an energy of 
200 MeV. 
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FIG. 3: The vector analyzing power An as in Fig. 2 except 
at very large scattering angles. 



MeV. At a scattering angle of six degrees An = —0.361 
ppm, see table I. This is comparable to A and thus a 
potential source of systematic error for the parity mea- 
surement. However, this error can be controlled by a 
measurement of An during the PREX experiment. 

Figure 5 illustrates how An depends on nuclear charge 
Z. It shows resuhs for ^^O, ^OCa, ^°Zi: and ^°^Ph at 850 
MeV. At forward angles the magnitude of An increases 
rapidly with increasing target Z. For ^^O, An is small 
except in the diffraction minima where the denominator 
in Eq. is small. For light nuclei the diffraction min- 
ima are sharp and deep. They become broader and less 
deep in heavy nuclei, see Fig. 1. This is reflected in the 
structure observed for An- 

Finally, HAPPEX He is measuring the parity violat- 
ing asymmetry A for elastic electron scattering from ^He 
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FIG. 4: The vector analyzing power A„ for elastic electron 
scattering from '^"^Pb versus scattering angle at an energy of 
850 MeV. 
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FIG. 5: (Color online) The vector analyzing power An for 
elastic electron scattering from ^'^O (blue dotted line), '*°Ca 
(green dashed line), ^"Zr (red dot-dashed line), and ■^°*Pb 
(black solid line) versus scattering angle at an energy of 850 
MeV. 



at 3 GeV and angles from six to nine degrees. To con- 
trol systematic errors, An will also be measured during 
HAPPEX He. Figure 6 shows our predictions for "^He at 
3 GeV using the three parameter Fermi charge density 
from ref. [T^. We find An to be very small, of order 
lO^^*^ and it passes through zero near 7.7 degrees. These 
results are so small compared to A because of the small 
Z of ''He, the high energy, and cancellations near the 
zero crossing. For completeness, we have collected in Ta- 
bles I and II our An predictions for the kinematics of the 
HAPPEX He and PREX experiments. 



TABLE I; Differential cross section and analyzing power A 
for elastic electron scattering from ^"^Pb at 850 MeV. 



e (deg. 


) da/dQ, (fm^) 






4 


1.76E-H03 


-1.14E- 


07 


4.2 


1.33E-H03 


-1.31E- 


07 


4.4 


1.00E-H03 


-1.49E- 


07 


4.6 


7.59E-I-02 


-1.69E- 


07 


4.8 


5.76E-I-02 


-1.90E- 


07 


5 


4.37E-f02 


-2.13E- 


07 


5.2 


3.33E-f02 


-2.38E- 


07 


5.4 


2.53E-f02 


-2.65E- 


07 


5.6 


1.93E-h02 


-2.94E- 


07 


5.8 


1.47E-h02 


-3.26E- 


07 


6 


1.12E-h02 


-3.61E- 


07 


6.2 


8.50E-h01 


-3.99E- 


07 


6.4 


6.48E-f01 


-4.41E- 


07 


6.6 


4.94E-I-01 


-4.85E- 


07 


6.8 


3.77E+01 


-5.33E- 


07 


7 


2.89E-h01 


-5.85E- 


07 


7.2 


2.22E-h01 


-6.39E- 


07 


7.4 


1.71E+01 


-6.94E- 


07 


7.6 


1.33E+01 


-7.50E- 


07 


7.8 


1.04E+01 


-8.05E- 


07 


8 


8.17E-H00 


-8.56E- 


07 


8.2 


6.51E-h00 


-9.01E- 


07 


8.4 


5.25E+00 


-9.37E- 


07 


8.6 


4.29E-I-00 


-9.64E- 


07 


8.8 


3.55E+00 


-9.79E- 


07 


9 


2.97E-H00 


-9.84E- 


07 


9.2 


2.51E-h00 


-9.78E- 


07 


9.4 


2.15E-(-00 


-9.64E- 


07 


9.6 


1.85E-(-00 


-9.44E- 


07 


9.8 


1.60E-h00 


-9.20E- 


07 


10 


1.39E-H00 


-8.93E- 


07 



TABLE 11: Differential cross section and analyzing power A 
for elastic electron scattering from '*He at 3 GeV. 



e (deg. 


) da/dQ. (fm^) 




5.88 


1.72E-02 


-1.29E-10 


6.08 


1.41E-02 


-1.25E-10 


6.28 


1.16E-02 


-1.18E-10 


6.49 


9.57E-03 


-1.09E-10 


6.69 


7.91E-03 


-9.81E-11 


6.9 


6.55E-03 


-8.41E-11 


7.1 


5.43E-03 


-6.72E-11 


7.31 


4.51E-03 


-4.73E-11 


7.51 


3.75E-03 


-2.40E-11 


7.71 


3.12E-03 


2.84E-12 


7.92 


2.60E-03 


3.35E-11 


8.12 


2.17E-03 


6.82E-11 


8.33 


1.81E-03 


1.07E-10 


8.53 


1.51E-03 


1.51E-10 


8.73 


1.26E-03 


2.00E-10 


8.94 


1.05E-03 


2.54E-10 


9.14 


8.80E-04 


3.13E-10 



6 




(degrees) 



FIG. 6: (Color online) The vector analyzing power An for 
elastic electron scattering from *He (solid line) versus scat- 
tering angle at an energy of 3 GeV. Also shown is A„ divided 
by 10 (red dashed line) to illustrate structure in the diffraction 
minimum. 



IV. CONCLUSIONS 

The analyzing power A„ gives the fractional difference 
in cross section for scattering to the left versus right for 
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